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基于 优先 队列 的 时 变 网 络 最 短路 径 算 法 
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摘 要 : 提出 了 基于 优先 队列 的 时 变 网 络 最 短路 径 算 法 ， 能 克服 传统 最 短路 径 算 法 难以 对 时 变 网 络 求解 最 短路 径 的 缺 
陷 。 提 出 的 时 间 窗 选择 策略 能 够 在 算法 求解 过 程 中 为 节点 选择 合适 的 时 间 窗 以 降低 路 径 长 度 ， 从 而 求 得 精确 解 。 进 一 
步 地 ， 算 法 使 用 了 优先 队列 组 织 节点 集合 以 提高 计算 效率 。 在 随机 生成 的 网 络 数据 以 及 美国 道路 数据 上 的 实验 表明 ， 
基于 优先 队列 的 时 变 网 络 最 短路 径 算 法 与 经 典 方法 相 比 ， 不 仅 能 够 来 得 精确 解 ， 运 算 速度 也 有 所 提高 。 
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Time varying network shortest path algorithm based on priority queue 
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Abstract: This paper presented a time-varying network shortest path algorithm based on priority queue to solve the time-varying 
network shortest path problem which was difficult for traditional shortest path algorithm. Proposed algorithm could address 
optimal solution by using proposed time-window selection strategy which could select appropriate time window for the node to 
reduce the path length. Also, the algorithm used the priority queue to organize node set, which could improves the computational 
efficiency. Experiments on randomly generated network data and united state road data show that compared with classical 
algorithm, the proposed algorithm is not only able to obtain global optimal solutions, but also improve the speed of the algorithm. 
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Throup00 使 用 优先 队列 解决 确定 网 络 下 的 单 源 最 短路 径 问题 ， 
利用 优先 队列 的 高 效 性 有 效 提 高 了 最 短路 径 求 解 的 速度 。 

最 短路 径 问题 是 图 论 和 算法 设计 中 的 一 个 基本 问题 ， 旨 在 千 统 算法 难以 解决 时 变 网 络 最 短路 径 问 题 。 时 变 网 络 与 普 
求解 给 定 网 络 中 给 定 两 点 之 间 的 最 短路 径 ， 在 路 径 规划 1、 通 的 静态 网 络 不 同 之 处 在 于 ， 一 点 到 另 一 点 的 距离 随 着 出 发 时 
GPS 导航 BY、 社交 网 络 5.l 等 领域 应 用 非常 广泛 。 间 的 变化 而 变化 。 这 也 就 使 得 时 变 网 络 中 的 最 短路 径 问 题 不 仅 
在 众多 的 最 短路 径 算法 中 ，Dijkstral 是 求解 单 源 最 短路 径 。 需要 考虑 路 径 中 的 节点 顺序 问题 ， 还 需要 考虑 节点 的 出 发 时 间 
的 代表 算法 ， 其 使 用 贪心 策略 能 够 有 效 求 得 确定 网 络 中 最 短路 。 问题。 如 图 1 所 示 ，weight 为 走 过 此 边 所 需 的 时 间 ，t 为 对 应 
径 。 但 由 于 时 间 复 杂 度 较 高 , Dijkstra 算法 并 不 能 在 大 规模 网 络 weight 有 效 的 时 间 区 间 。 若 求 从 节点 1 至 节点 3 的 最 短路 径 ， 
中 表现 出 良好 的 性 能 。 为 解决 这 一 问题 ， 学 术 界 对 其 进行 了 大  ” 则 有 可 选 路 径 1 一 2 一 3 和 1 一 3， 且 两 条 路 径 上 分 别 有 4 个 和 2 
量 研究 。 王 华中 利 用 邻接 点 算法 对 Dijkstra 算法 进行 改进 , 在 个 时 间 窗 。 如 何在 算法 运行 过 程 中 动态 确定 路 径 节点 ， 并 正确 
定 程度 上 降低 了 算法 的 时 间 复 杂 度 和 空间 复杂 度 ; 马 小 雨 等 人 ”选择 时 间 窗 以 使 路 径 最 短 ， 是 传统 方法 求解 时 变 网 络 最 短路 径 
中 提出 从 减少 搜索 计算 顶点 数量 入 手 ,对 经 典 Dijkstra 算法 进行 ”的 原因 。 
改进 ， 提 高 了 算法 运行 的 效率 ; 王 树 西 等 人 09 针 对 算法 无 法 对 时 变 网 络 比 普通 网 络 更 为 符合 生产 生活 实际 ， 在 公路 导 骨 
不 连通 有 向 图 等 问题 进行 了 改进 ,并 将 其 用 于 道路 选择 系统 中 。 023、 卫星 网 络 03 和 路 径 规划 0 等 方面 都 更 符合 实际 情况 。 例 如 ， 
尽管 上 述 改进 算法 具有 很 多 优点 ,但 其 仍 存在 着 许多 约束 条 件 ， ”在 实际 的 道路 交通 环境 下 ， 一 条 道路 在 不 同 的 时 间 区 间 下 的 通 
未 能 解决 高 响应 、 大 规模 环境 下 的 网 络 最 短路 径 问 题 .2004 年 ， ”过 时 间 可 能 是 变化 的 ， 这 时 ， 使 用 时 变 网 络 模型 来 模拟 道路 交 
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通 更 加 科学 。 学 术 界 已 对 动态 网 络 最 短路 径 问 题 进 行 了 广泛 研 
究 ， 并 获得 了 一 定 的 研究 成 果 。 戴 浴 琴 等 人 (3 使 用 改进 的 
Dijkstra 算法 解决 卫星 时 变 网 络 中 的 最 短路 径 算 法 , 邹 亮 等 人 0 
使 用 遗传 算法 解决 动态 网 络 中 的 最 短路 径 问 题 ， Sever 等 人 09 
使 用 以 动态 规划 为 基础 的 综合 策略 的 方式 解决 变动 网 络 中 的 最 
短路 径 问 题 。 以 上 方法 中 大 多 需要 进行 多 次 遍历 ， 因 此 存在 效 
率 低 下 的 问题 ,基于 启发 式 算法 还 存在 求解 精确 度 不 足 等 问题 。 


te[0,6) weight =5 

tel[6,10]weight =4 
te[0,3)weight = 4 
tel[3,5]weight =5 


ft e[0,5)weight =5 
tel[6,9]weight = 2 


图 1 时 变 网 络 例 


本 文 提出 了 一 种 基于 优先 队列 的 时 变 网 络 最 短路 径 算 法 ， 
法 使 用 时 间 窗 选择 策略 ， 同 时 运用 优先 队列 高 效 性 的 优点 ， 
够 快速 求解 时 变 网 络 最 短路 径 问 题 。 在 随机 时 变 网 络 数据 与 
国道 路 地 图 上 的 实验 表明 ， 提 出 的 算法 与 经 典 算 法 相 比 ， 不 
能 求 得 精确 最 短路 径 ， 求 解 速度 也 有 所 提高 。 
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1 ”预备 知识 


优先 队列 53 是 一 种 基于 堆 的 动态 排序 集合 ， 
又 堆 为 基础 进行 构造 。 利 用 其 在 数据 结构 上 的 优势 ， 能 够 有 效 
组 织 节 点 并 提高 计算 速度 。 优 先 队列 五 支持 如 下 操作 : 

findMin(H) : 返回 H 中 最 小 的 元 素 。 由 于 在 优先 队列 中 最 
小 元 素 在 队列 首位 ， 返 回 队 列 中 第 一 个 元 素 即 可 。 此 操作 时 间 
复杂 度 为 0() 。 

insert(H,q) :向 右 中 插入 元 素 4。 由 于 在 插入 后 需要 对 优先 
队列 进行 维护 ,此 操作 时 间 复 杂 度 为 cdogm) 。 

dec 一 key(H,a,*) :将 对 象 4 的 值 降 至 x*。 若 a 的 值 小 于 *, 则 返 
可 一 个 错误 。 此 操作 时 间 复 杂 度 为 odognn)。 

extractMin(H) :去 掉 并 返回 # 中 最 小 值 的 元 素 , 即 队列 首 元 
素 ,此 操作 时 间 复 杂 度 为 ologn) 。 

优先 队列 主要 应 用 于 基于 事件 驱动 的 模型 中 ， 一 般 用 于 对 
操作 系统 中 进程 的 管理 。 操 作 系 统 将 进程 按 优 先 级 或 执行 时 间 
进行 排序 后 ， 可 每 次 从 队列 首位 取得 需要 操作 的 进程 。 得 益 于 
其 较 高 的 排序 效率 ， 在 队列 中 某 个 进程 的 优先 级 发 生变 化 时 ， 
先 队列 能 快速 地 对 队列 进行 重新 排序 。 


常 使 用 最 小 二 


进行 处 理 ， 即 确定 出 发 时 间 ;， 在 改变 其 他 节点 与 开始 节点 的 距 
离 时 ， 快 速 地 重新 排序 。 因 此 ， 利 用 优先 队列 进行 时 变 网 络 最 
短路 径 规 划 是 可 行 的 。 


2 问题 定义 


定义 1 时 变 网 络 (time varying network, TVN) 
给 定 一 个 有 向 连通 图 G(V,A,W) , 其 中 了 = fw) 为 图 中 


节点 集合 ，4={Q(,v))]v,vjeV} 为 图 G 中 边 的 集合 ， 丈 为 一 组 时 
变 函 数 ， 则 称 G 为 时 变 网 络 。 

在 时 变 函 数 W(DeW 中 ，W(D) 为 自 记 至 六 且 出 发 时 刻 为 
1 时 ， 边 www) 的 长 度 。 

定义 2 ”时间 窗 (time window) 

在 时 变 网 络 中 , 时 间 区 间 [,) 称 为 时 间 窗 。 其 中 1 为 时 间 窗 
下 界 , w 为 时 间 窗 上 界 。 对 于 时 间 窗 集合 区 ,ma [zw) ,车 1=0 
且 A<k = (=D.…m) ， 则 称 此 集合 为 时 间 窗 集 。 

设 边 ww) 上 有 时 间 窗 集 mw ，mm 包含 严 个 时 间 窗 ， 即 


ws oo |i=1.… 册 , 则 在 时 变 网 络 中 ， 有 


0 tem, 
W,,(D) =4... (1) 
0, fetw, 
其 中 : 9 为 时 间 1 属 于 在 时 间 窗 tw 中 时 ， 边 ww) 的 长 度 。 
定义 3 ”时 变 网 络 的 路 径 (a path of TVN) 


对 于 有 序 序列 path=< G4),…,v,,4,)> ,vw 为 序列 中 第 i 个 节 
点 ,5 为 vi 的 出 发 时 间 ， 则 称 path 为 时 变 网 络 的 一 条 路 径 。 由 
此 也 可 以 得 出 ， 路 径 path 的 节点 路 由 为 节点 序列 <wm mm > 。 

时 变 网 络 最 短路 径 求解 问题 即 在 给 定 的 时 变 网 络 G 中 为 
每 个 节点 v 选择 全 局 最 优 的 出 发 时 间 : ,使 自 起 点 8 至 终点 7 的 


距离 最 小 。 由 此 ， 时 变 网 络 的 最 短路 径 问题 可 定义 为 
min 2 Xpa * 4 二 ) 
subject to: 区 2 一 六 xi =1 
(LpP)sA (gD)eA 
2 Xpg > Xap =0 (2) 
(p.g)eA (gq,p)eA 


x,, E(0,D,v(p,g)eA 


pq 


其 中 ，> XxxW(t) 为 最 短路 径 长 度 ，*j 为 边 (p,q) 是 否 存 


在 于 网 络 中 的 标志 ; 


Wu 为 边 (p,9) 在 5 时刻 出 发 的 长 度 。 
3 ”基于 优先 队列 的 时 变 网 络 最 短路 径 算法 
在 时 变 网 络 中 ， 两 点 之 间 的 距离 由 起 点 的 出 发 时 间 决 定 。 


在 时 变 网 络 最 短路 径 规划 问题 中 ， 由 于 两 点 之 间 的 长 度 


因此 要 求解 时 变 网 络 中 的 最 短路 径 ， 应 逐个 确定 节点 的 出 发 时 


前 一 个 节点 的 出 发 时 间 决 定 ， 且 一 条 边 长 的 确定 会 影响 周围 多 
条 边 的 长 度 ， 所 以 对 于 时 变 网 络 最 短路 径 规 划 问 题 ， 需 要 逐个 


间 。 所 以 ， 算 法 主要 问题 在 于 如 何 高 效 组 织 节 点 ， 并 在 对 迭代 
过 程 中 节点 的 变化 及 时 进行 维护 。 优 先 队列 具有 便于 取出 元 素 、 


确定 节点 的 出 发 时 间 ， 并 在 确定 之 后 更 新 相 邻 边 的 长 使 用 


VE 


排序 快 、 维 护 代价 低 等 优点 。 采 用 优先 队列 按 与 出 发 节点 的 距 


tC 先 队列 管理 网 络 节点 的 方法 和 效果 与 进程 管理 相似 首先 按 
节点 距离 路 径 开始 节点 的 距离 进行 排序 ， 每 次 取 队 列 首位 节点 


离 的 方式 组 织 节 点 ， 便 于 在 途 代 过 程 中 快速 取得 与 出 发 节点 最 
近 的 节点 。 在 更 新 过 程 中 ， 若 某 节点 的 情况 发 生变 化 ， 优 先 队 
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列 可 快速 调 
织 时 变 网 络 
3.1 算法 流 


整 其 位 置 ， 维 护 代价 较 低 。 因 此 ， 采 | 
中 的 节点 更 为 高 效 。 


i 程 


图 2 给 


的 时 变 网 络 


合 出 了 算法 的 流程 。 


算法 基本 流程 如 下 : 


进行 处 理 
的 后 置 节 点 选择 时 间 窗 并 更 新 其 与 出 发 节点 的 


队列 首 节 点 
ed 
开始 ， 沿 前 
法 可 分 为 三 
生成 最 短路 
QO 优先 
优先 队列 ， 


列 首 节点 ， 直 至 移 除 了 目标 节点 ; 最 后 从 
点 向 出 发 节点 回 淹 

ey 即 优先 队列 初始 化 部 分 、 
径 部 分 。 算 法 的 各 个 部 
队列 初始 化 。 负 责 对 给 


分 功能 分 述 如 下 : 


目标 节 
， 得 到 最 短路 径 。 据 此 ， 算 
节点 更 新 部 分 及 


优先 队列 组 


首先 对 给 定 
， 得 到 优先 队列 ; 然后 进行 迭代 循环 ， 为 
距离 ， 
让 点 


全 定时 变 网 络 进行 搜索 以 生成 


初始 化 节点 。 在 节点 搜索 过 程 中 ， 首 先 使 用 


多 先 搜索 算 
优先 队列 及 
节点 5 作为 
初始 化 。 

回 节 点 


法 获得 有 序 节点 
节点 初始 化 。 在 
搜索 的 起 点 ， 以 


序列 ， 
广度 搜索 过 程 中 ， 使 用 


给 定 的 


更 新 。 负 责 节 点 返 代 更 新 及 优先 队列 维 


最 短路 径 长 度 


间 窗 选择 及 
重复 此 过 程 


LE 护 ， 
。 每 次 迭代 首先 对 队列 首 节点 的 后 
相关 信息 更 新 ,然后 将 队列 首 节点 移 除 并 维 
， 直 至 将 目标 节点 自 优先 队列 移出 。 


加 生成 


最 短路 径 。 


录 当 前 最 短 


路 径 中 本 节点 的 前 置 节点 


并 为 每 个 节点 赋 初 值 ， 完 


以 求 日 
节点 进行 时 
护 队 列 。 


使 S 为 队列 首 节点 。 最 后 对 5 进行 


上 


LT 


在 节点 更 新 的 过 程 中 ， 会 同时 为 节点 记 
。 在 节点 更 新 完成 后 ， 通 


自 目标 节 
路 径 。 


过 


QD 优先 队 
列 初始 化 


© 


节点 更 
新 迭代 


点 的 前 置 节点 向 出 发 节点 递归 下 


优先 队列 及 节点 初始 化 


初始 化 出 发 节 


为 队列 首 节 点 的 后 置 节 点 选 
择 时 间 窗 
vy 
更 新 后 置 节 点 相关 值 

维护 优先 队列 


vy 
移 除 队列 首 结 点 
维护 优先 队 罗 


并 


标 节点 向 出 发 节 ; 点 可 济 
得 到 最 短路 径 


y 
( 结束 ) 


图 2 基于 优先 队列 的 时 变 网 络 最 短路 径 算 法 流程 


3.2 算法 形式 化 描述 


介绍 算法 的 形式 化 接 


中 的 符号 及 


。 在 开始 


符号 释义 ， 见 表 1。 


算法 


滴 ， 即 可 得 到 最 短 
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表 1 符号 描述 

符号 符号 含义 
G 时 变 网 络 
A 符合 优先 队列 结构 的 节点 集合 
BS(v) 节点 v 的 后 置 节 点 
L 己 选择 时 间 窗 的 节点 集合 
S 出 发 节点 ， 即 最 短路 径 起 点 
了 标 节 点 ， 即 最 短路 径 终点 
tw 边 (p,q) 上 的 时 间 窗 集 
Pw tw 的 第 个 时 间 窗 
1 twi, 的 下 界 
uy wi 的 上 界 
0 选择 tw 时 边 (p,q) 的 长 度 
下 vy 的 父 节点 
t, 7 出 发 时 间 
d, 点 v 当前 最 短路 径 长 度 ， 即 与 8 的 距离 


于 对 应 关系 如 下 : 步 又 


基于 优先 队列 的 时 变 网 络 最 短路 径 算法 如 下 所 示 。 其 与 流 
1 一 3 为 优先 队列 初始 化 部 分 ， 步 又 


4 一 6 为 节点 更 新 部 分 ， 步 骤 7 一 8 为 最 短路 径 生 成 部 分 。 
算法 1 基于 优先 队列 的 时 变 网 络 最 短路 径 算法 


输入 : G、S、T。 


令 


出 : 最 短路 径 path 及 最 短路 径 长 度 dy 
5 为 起 点 ， 对 图 G 进行 广度 搜索 ,获得 包含 图 中 全 部 节点 的 有 序 序列 
={5,.…,4,.…T,.…4,} ， 设 集合 L={@)} 
4 中 的 每 个 节点 v 设 d=%， 环 = 
dso =0,， 即 ds=0 
Tg 工时 ， 循 环 步 又 4-5: 
算法 2 更 新 4[0] 及 BS(A[0]) 
工 = {LL, A[0]} ， 并 执行 extractMin(A) 


执行 算法 3， 得 到 最 短路 径 path 


输 


痢 为 
后 对 


不 重 
因此 
除 无 
应 的 


为 此 
ee 
则 等 


首 节 上 


出 ， 得 到 最 短路 径 path 及 最 短路 径 长 度 di 

算法 2 为 节点 更 新 算法 。 此 子 算法 首先 使 用 时 间 窗 选择 策 
节点 选择 时 间 窗 ， 而 后 更 新 节点 4 的 当前 最 短路 径 4, ， 最 
亿 先 队列 进行 维护 。 在 算法 输入 部 分 中 ， 4[0] 为 优先 队列 
点 ， BS(4[0] 为 时 变 网 络 中 4[0] 的 后 置 节 
第 2 意 定 义 2 及 3 可 知 ， 时 间 窗 集合 是 由 一 系列 非 负 且 
县 的 时 间 区 间 组 成 的 ， 每 个 时 间 区 间 都 有 其 对 应 的 边 长 。 
在 为 节点 选择 时 间 窗 时 ， 首 先 根据 前 置 节点 的 出 发 时 间 排 
效 的 时 间 窗 ， 然 后 计算 每 个 时 间 窗 在 其 最 早 时 间 出 发 的 对 
路 径 长 度 ,最 后 在 路 径 长 度 最 小 的 时 间 窗 中 取 其 出 发 时 间 。 
根据 分 析 ， 可 得 出 时 间 窗 选择 策略 : 
为 mw 中 的 每 个 时 间 窗 计算 对 应 长 
时 间 窗 过 期 ， 令 d)=%; 
出 发 令 d)=d,+0,; 


点 


度 ad : 若 d, >w ， 则 视 
若 d, em ， 即 4, 在 时 间 窗 中 ， 
4d, ， 即 d, 早 于 此 时 间 
"= ,+0, 


pq 9 


车 4 < 窗 ， 
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3.3 ”算法 实例 


节点 7 


录用 稿 
为 使 得 节点 4 能 取得 最 短 长 度 , 车 4d, <d,，, 则 使 4,=4d,，, 并 


记录 当前 最 短路 径 的 前 置 节 点 f=P 。 


算法 2 节点 更 新 算法 


输入 : AI0] 、BS(A[0D、twj、d,、d,。 
输出 :更 新 后 的 d 及 更 新 后 的 优先 队列 4 。 
令 p=AI0] ，qe BS(A[0]) 对 于 npw 中 的 每 个 时 间 窗 tw : 
若 d, < ， 则 设 忆 = Os te 


车 lb <d < upg i 则 设 d=d, +0% t,=d 


p 


车 4d, >wW，。 ”， 则 设 d) =% 


车 d,>d， ， 则 设 d=4d),， =p ,t= 


对 A 执行 dec 一 key(A,g,d,) 
算法 结束 
当 节 点 迭代 完成 后 , 集合 ZL 中 除 5 以 外 的 所 有 节点 的 五 都 
己 被 赋值 。 可 使 用 算法 3 从 7 开始 向 前 根据 到 向 前 迭代 回溯 至 
3 得 到 所 求 最 短路 径 。 
算法 3 ”路径 回 溯 算 法 
输入 T、f,、S。 


出 : 最 短路 径 path 。 
设 g4=7 ， 并 有 有 序 序列 path =< gq > 
循环 步骤 3-4 


当 5S gpath 时 ， 


path =< gq, path > 
算法 结束 ， 得 到 最 短路 径路 


path 


为 了 进一步 说 明 算法 的 运行 过 程 ， 下 面 使 用 图 3 中 的 时 变 
网 络 进行 算 例 演示 。 图 3 中 共有 5 个 节点 ， 出 发 节点 5 和 目标 


已 经 标 出 。 


时 变 网 络 的 详细 信息 见 表 2。 


网 络 中 的 每 条 边 都 有 1 一 2 个 


时 间 窗 。 其 中 时 间 窗 为 [0,%) 的 表示 不 论 何 时 出 发 , 边 的 长 度 都 
为 对 应 的 9 值 。 
图 4 将 算法 运行 过 程 以 优先 队列 结构 进行 展示 ， 圆 内 标号 


为 节点 序号 , 插 号 内 为 节点 4,。 图 4 中 四 为 初始 化 后 的 优先 队 
列 结构 ， 从 @ 开 始 使 用 算法 2 对 优先 队列 进行 更 新 ，@ 为 最 短 
路 径 长 度 及 使 | 


j 算 法 3 求 得 的 最 短路 径 。 
F 面 对 每 个 步骤 图 进行 解释 ; 
运行 算法 1 的 第 1 一 2 步 。 使 用 ) 


度 搜 索 对 图 G 进行 搜索 ， 
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生成 优先 队列 4 ， 
最 后 初始 化 5 点 。 


刊 


Chi 
Re 


对 4 中 的 节点 进行 初始 化 得 到 优先 队列 ， 


表 2 时 变 网 络 信息 

节点 后 置 边 1 0 
[0,D 1 

(S,1) 
S [Loo) 3 
(S,2) [0,2) 1 
1 (1,T) [0,%) 5 
(23) [0,%) 3 
2 [0,2) 4 

(2,T) 
[2,%) 1 
3 (3,T) [0,oo) 2 


对 4[0] 即 S$ 点 的 邻居 节点 1、2 使 


算法 2 更 新 节点 信息 


其 中 时 间 窗 ws 及 ws; 都 选择 了 第 1 个 时 间 窗 ，2 均 为 1。 


将 队列 首 元 素 8 并 入 工 中 ,并 从 4 中 删除 ,维护 优先 队列 。 
此 时 队列 首 节点 为 节点 1。 
运行 算法 2 更 新 节点 1 的 邻居 节点 7 并 维护 优先 队列 。 更 


St 


后 必 =6 ， 由 于 dr < ， 


7 节点 在 队列 中 的 位 置 高 于 节点 3。 


将 队列 首 元 素 节 点 工 并 入 工 中 ， 从 4 中 


| 除 节点 1 并 维护 


优先 队列 。 此 时 队列 首 节点 为 节点 2。 

运行 算法 2 更 新 节点 2 的 邻居 节点 3、7， 
了 wi =[2,%), di=3, d=4. 
将 队列 首 元 素 节 点 2 并 入 L 中， 并 从 4 中 
队列 首 节点 。 


更 新 后 7 


元 素 7 并 入 工 中 ， 并 从 优先 队列 中 
束 ， 得 到 最 短路 径 长 度 必 =3。 


! © 1 加 | 国 | 
1 -一 上 1 二 

(s ‘0)) ' Go) (10) | 
| Co | ~ Nt | > ~ | 
| Ne) ( 0 ro ) (0 1 (sg ) 20) 
1 SA ed = 和 | 
1 es ~ I 
Ce mm ‘0) I 
/A 人/ NN7/ RD J | 
人 
1 5 1 上 | 
I | | 
1 \ 7/ 1 J 1 -一 | 
1 -天 SO | (20 ) | (am) | 
| 79 ) E09) | | > | 
1 Ne ee 1 1 一 一 | 
| | (re ) 9) | (10 ) | 3 ) | 
a) | A A NA WN ' 

1 1 
| 一 L={S(0)} | L={S(0),1(D} | L={S(0),1D} | 
1 '@®@ 1'@ | 
Ar 一 1 _ 1 

| i 上 pi 1 | 
| Wy ( 300 ) d; =3 
| py se | Dy | path= | 
| 9) | | {5,0),(2,D,07,3)} | 
1 本 1 1 | 
| L={S(0),1(1),2(1)} | L={5(0),1D,20D TO)}! | 
1 1 1 | 
下 二 OE ER | 


图 4 算 例 计算 过 各 


由 于 7 没有 后 置 节 点 ， 运 行 算 法 2 后 无 值 更 新 。 将 队列 首 
出 除 。 此 时 节点 循环 更 


选择 


删除 。 此 时 了 为 


新 结 


录用 稿 
使 用 算法 3 完成 最 短路 径 回 滴 ， 得 到 最 短路 径 path ， 
结束 。 


表 3 给 出 了 算法 运行 过 程 中 每 个 节点 各 项 值 的 变化 。 关 注 
标 节点 7 在 图 4 的 由 和 人 @@ 的 两 次 数据 变化 ， 可 以 发 现时 间 窗 


选择 策略 在 寻找 最 优 解 中 发 挥 的 作用 : 
在 由 中 ， 节点 1 为 队列 首 元 素 ， di=1 o 


间 窗 四 =[0,oo) Or =5 ， 


时 间 窗 wl; =[0,2) 2 =4 和 m3 =[2,%) 03, =1 


因此 a; =d +O7 =0 。 
在 @ 中 ， 节 点 2 为 队列 首 元 素 ，d,=1。 
, 对 两 个 时 间 窗 分 析 


算法 


边 47) 上 有 唯一 时 


边 (27 上 有 2 个 


若 选择 ow; ， 即 在 1 时 刻 立 即 出 发 ， dv = 玫 + 同 =1+4-5， 
若 选择 mx ， 即 在 2 时 刻 出 发 ，dyr = 尽 +Or =2+1=3。 


所 以 ， 选 择 mzr ， 即 在 节点 2 处 等 
时 刻 2 出 发 ， 可 以 使 得 di 达到 最 小 值 3。 


表 3 节点 状态 


待 1 个 单位 时 间 ， 并 在 


图 号 节点 4 下 tm 所 在 集合 
S 0 人 人 A 
1 o 人 人 4 
@ 2 o 人 亿 A 
3 o 人 人 4 
了 o 人 人 A 
S 0 人 人 L 
1 1 S [0.1) A 
©-® 2 1 S 0,2) A 
3 o 人 人 A 
了 o 人 人 4 
S 0 @ 人 L 
1 1 S [0,D 元 
@-© 2 1 S [0,2) A 
3 o 人 人 A 
T 6 1 [0,%) A 
S 0 人 人 L 
1 1 S [0,D L 
©-O 2 1 S [0,2) L 
3 4 昂 [0,%) A 
了 3 多 [2,%) A 
S 0 人 人 L 
1 1 S [0,D 万 
® 2 1 S [0,2) L 
3 4 2 0,%) A 
注 : 图 中 深 色 格 为 迭代 后 变化 的 内 容 


3.4 算法 分 析 
3.4.1 正确 性 分 析 
下 面 使 用 反 证 法 来 证 明 算 法 的 正确 性 。 


法 实例 可 以 发 
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现 ， 在 算法 运行 过 程 中 ， 任 意 时 刻 时 变 网 络 中 的 任意 一 个 节点 
i 合 合 L。 由 此 有 以 下 定理 : 
上 的 所 有 节点 都 属于 集合 L 。 

证 明 由 图 5 所 示 ， 假设 有 从 5 到 7 的 路 径 了 和 路 径 p'， 
其 中 p'=pi+p,，, 其 中 P 中 的 节点 全 部 属于 集合 L ,在 p' 中 有 至 
少 一 个 节点 在 集合 4 中 ， 即 图 中 的 a 点 。 当 算法 停止 时 ，T eL 
而 ae 工 ， 说 明 有 dr<4d。。 同 时， 由 于 图 中 边 长 皆 为 正 值 ， 所 以 
有 ?的 长 度 小 于 p'。 

里 1 得 证 。 

in, 

面 定理 

策略 能 为 节点 选择 最 优 时 间 


度 证 明了 最 短路 径 中 的 节点 均 已 求 
2 从 时 间 窗 角度 证 明了 时 间 窗 选择 
窗 以 求 得 当前 最 短路 径 。 


LI1A 


图 5 正确 性 说 明 图 


定理 2 最 短路 径 上 的 任 
路 径 长 度 都 是 最 小 的 。 

证 明 设 有 最 短路 径 path=< (4),…,(v,h)> ， 对 于 任意 节 
点 ,其 当前 路 径 长 度 di; 于 0 由 决定 , 因此 di 主 
要 受 t 影响 。 下 面 讨论 二 在 不 同情 况 下 对 di; 的 影响 : 

(1) 设 有 <t ,使 得 d; 算法 2 的 步骤 1 可 知 ， 
二 可 取 值 可 能 为 & 或 六 〈 选 定时 间 窗 下 界 )。 若 # <d ,， 根 


意 节 点 的 出 发 时 间 所 对 应 的 当前 


=1,+0 ,I 


=t+0 :1 


据 算法 2 步骤 1.c，du =%; 若 <liin ,根据 算法 2 步骤 1.d， 


d; <d; 。 因 此 ， 若 ;<4 则 定 有 d; <d;。 
(2) 设 有 4 >t， 使 得 d;=#+0G， 则 
与 {在 同 
间 窗 内 , 则 有 di 
车 #6 与 4 不 在 同 


有 以 下 两 种 情况 : 
时 间 窗 内 , 或 不 在 同一 时 间 窗 内 。 车 与 在 同一 时 
f;+E+0 ,由 于 s 为 任意 正 数 , 所 以 d; <d,; 
时 间 窗 内 ， 根 据 算 法 2 步骤 1.4d，d; <d;。 因 


此 ， 若 > 则 定 有 di<di。 
定理 2 得 证 。 
3.4.2 时 间 复 杂 度 分 析 


定理 3 在 具有 7 个 节点 和 区 个 时 间 窗 的 时 变 网 络 中 ，] 
于 优先 队列 的 时 变 网 络 最 短路 径 算 法 的 时 间 复 杂 度 为 
o(3n+mn+nlog(n)) 。 


证 明 下 面 对 三 个 算法 


及 


进行 时 间 复 杂 度 分 析 : 
算法 2: 假设 每 个 节点 有 太 个 时 间 窗 ， 步 又 1 需要 循环 尹 
次 ， 步 又 2 对 优先 队列 的 维护 在 最 差 情况 下 需要 odog(n)) 的 时 
间 。 所 以 算法 2 的 总 时 间 为 m+o(log(n)) 。 
算法 3: 此 算法 是 一 个 递归 算法 ， 用 于 递归 求 出 


最 短路 径 
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录用 稿 


路 
间 复 杂 度 为 o(n) 。 


算法 1 对 节点 的 广度 搜索 及 节点 初始 化 需要 遍历 所 
点 ,即时 间 复杂 度 为 om)x2 ; 初始 化 出 发 节点 时 间 复杂 度 为 00D) ; 


步骤 4-6 的 循环 在 最 差 情 况 下 需要 循 


此 可 得 出 算法 的 时 间 复 杂 度 为 


o(n)+o(n)+o() +nx(m+o(log(n))) +o(n) 


=0(3n+n’ +nlog(n)) 。 


定理 3 得 证 。 


定理 3 可 得 到 如 下 关于 时 间 


若 m>n ， 算 法 的 时 间 复 杂 


/ 
显影 响 ; 


数量 与 算法 性 能 的 推论 : 
推论 1 若 m<log(n) ,时 间 窗 数量 将 不 会 对 算法 性 能 造 
度 将 显著 提高 。 


证 明 由 算法 的 三 个 主要 部 分 分 别 是 初始 化 、 时 间 窗 选择 、 


节点 迭代 ， 分 别 耗 时 oC2n) 、olmn) 和 o 


于 选择 时 间 窗 的 时 间 。 


(2) 车 mm>log(m) ， 则 olmm) >olnlog(n)) ,时 间 窗 选择 的 时 


将 超过 迭代 操作 的 时 间 。 


(1) 若 m<log(m) ， 则 olmn)<olnlog(n)) ， 算 法 迭代 的 时 


(nlog(n)) 。 


。 在 最 坏 情况 下 ， 所 有 节点 都 在 最 短路 径 上 ， 算 法 4 的 时 


有 贡 


环 寺 次， 所 以 时 间 复 杂 度 
为 nx(m+o(log(m)) ; 循环 结束 后 的 算法 4 时 间 复 杂 度 为 0(n) 。 


成 明 


间 大 


间 


(3) 车 m>n ， 时 间 复 杂 度 将 大 于 oG3n+r +nlog(m)) ， 时 间 


复杂 度 明显 


首长 。 此 时 ， 节 点 的 时 间 窗 复杂 度 会 较 大 地 影 


法 的 性 能 。 
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实验 


本 章 将 用 Dijkstra 算法 、Thorup [1 
的 时 变 网 络 最 短路 径 算法 在 5 个 随机 9 


1 的 方法 和 基于 优先 


响 算 


队列 
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表 5 随机 生成 地 图 下 的 实验 结果 
Dijkstra Thorup (2004) 提出 的 方法 
数据 编号 

Error% CT(ms) Error% CT(ms) Error% CT(ms) 

D， 10.74% 18 10.74% 7 0% 13 

D, 17.05% 136 17.05% 42 0% 46 

D;, 19.88% -377 19.88% 131 0% 138 

D， 23.90% 839 23.90% 226 0% 233 

D; 23.96% 1565 23.96% 366 0% 372 


表 5 为 三 种 算法 在 随机 网 


络 上 的 实验 结果 ， 包 括 每 种 算法 


在 不 同 数据 集 下 的 错误 率 以 及 运算 时 间 (CT)。 错 误 率 为 0 表 
示 算 法 达到 了 全 局 最 优 解 。 
< 用 了 优先 队列 结构 ， 在 5 组 


由 表 5 可 以 看 出 ， 提 出 的 方法 由 于 
数据 中 求解 时 间 增 长 幅度 不 大 ， 


林 
符合 时 间 复 杂 度 预 共 


; 相 比 之 下 , Dijkstra 算法 的 运算 时 


选择 过 程 ，Thourp 


让 一 
又 


; 数据 集 下 耗 时 是 


的 运算 速度 较 提 出 的 方法 稍 快 ; 
面 ， 提 出 的 算法 采用 
窗 并 求 得 最 短路 径 ， 在 5 组 实验 中 均 能 够 能 求 得 最 优 解 。 而 


Dijkstra 和 Throup 上 


于 直接 选用 


间 增 长 
提出 算法 的 4 倍 ; 由 于 没有 时 间 窗 
在 运算 精确 
首 ， 可 取得 最 优 时 间 


了 时 间 窗 选择 策 


首 个 可 用 时 间 窗 ， 从 而 产生 了 


一 定 且 相 同 的 误差 。 


4.2 美国 道路 地 图 测试 


为 了 测试 算法 在 真实 网 络 中 的 性 能 ， 


本 节 使 


用 EE 


国道 路 地 


E 成 的 时 变 网 络 及 美 


路 网 络 进行 实验 对 比 。 实验 环境 : Intel 
16 GB 内 存 ，128 GB SSD 硬盘 ， 
Studio 2015 编程 软件 ，C# 语 言 实现 。 
4.1 随机 地 图 测试 


本 组 实验 数据 共 5 组 ， 节 点 个 数 分 别 为 1000 一 9000 个 ,使 
] Randgraph (http://www.dis.uniromal.it/challenge9/download. 
shtml) 生成 ， 对 于 地 图 中 的 每 条 边 ， 参 考 此 边 的 长 度 随 机 添加 


RY Core™ i7-6700 3.4 


时 间 窗 。 表 4 概括 了 实验 网 络 结构 信息 。 


国道 
GHz， 


Windows 10 操作 系统 ，Visual 


图 (http:/www.dis.uniromal .it/challenge9/downloa d.shtml ) 进行 
测试 ,由 于 美国 道路 地 图 是 静态 道路 网 络 ,为 了 完善 实验 过 程 ， 
实验 前 为 网 络 中 的 部 分 边 添加 了 时 间 窗 。 表 6 为 实验 网 络 结构 
信息 。 
表 6 实验 国道 路 网 络 结构 信息 
实验 数据 节点 数量 边 数 量 时 间 窗 数量 
New York- D, 264,324 733,846 1,731,876 
Colorado- De 435,666 1,057,066 2,621,523 
Florida- Dr 1,070,376 2,712,798 6,917,634 


与 随机 网 络 数据 的 实验 方法 相同 ， 表 7 中 列 出 了 在 美国 i 
中 的 实验 结果 。 


路 地 图 


在 大 规模 网 络 中 ， 提 出 的 算法 依然 能 够 


求 得 全 局 最 


表 4 随机 网 络 结构 信息 
实验 数据 节点 数量 边 数 量 时 间 窗 数 和 
D， 1000 1507 4022 
D， 3000 4479 8743 
D， 5000 7450 16871 
D, 7000 10489 22354 
D; 9000 13542 30874 


每 组 实验 都 会 随机 选择 一 个 出 发 节点 和 一 个 目标 节点 


行 多 次 重复 实验 。 在 实验 过 程 中 ， 
短路 径 算 法 采用 时 间 管 
时 间 窗 以 决定 两 点 间 的 距离 。 


可 | 


基于 优先 队列 的 时 变 网 
选择 策略 ， 其 余 两 算法 均 采 | 


选择 


， 进 
络 最 
首 个 


亿 解 ， 而 其 他 算法 由 于 难以 选取 最 优 时 间 窗 ， 出 现 
了 10% 一 20% 的 计算 误差 。 计 算 时 间 上 ， 在 美国 道路 地 图 中 ， 


Dijkstra 算法 的 计算 时 间 上 升 最 快 ，26 万 节点 的 平均 计算 时 间 


为 9 min，100 万 节点 的 计算 时 间 
方法 在 26 万 节点 中 的 计算 时 间 在 3 s 以 


已 经 达到 了 约 1 h。 而 提出 的 
，100 万 节点 的 平均 


mm 


计算 时 间 约 为 9s， 体 现 出 优先 队列 结构 在 运算 速度 上 的 优势 ; 


时 间 窗 选择 会 消耗 一 部 分 运算 时 
Thorup0110% 的 时 间 换 取 了 11.37% 的 准确 度 ， 在 D; 组 中 用 


间 , 在 Dy 组 中 ， 


提出 的 算法 用 


Thorup003% 的 额外 运算 时 间 换 取 了 20% 的 准确 度 。 


结束 语 


本 文 提 出 了 一 种 基于 优先 队列 的 时 变 网 络 算法 用 


于 求解 时 
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录用 稿 


变 网 络 中 的 最 短路 径 问题 。 提 出 的 算法 采 / 


时 间 窗 选择 策略 以 


克服 传统 算法 难以 对 时 变 网 络 求 最 短路 径 精 确 解 的 问题 ,同时 ， 


算法 采 上 
机 实验 数据 及 真实 网 络 数据 的 实验 表明 ， 


优先 队列 组 织 网 络 节 点 ， 提 高 了 算法 求解 速度 。 在 随 
基于 优先 队列 的 最 短 


路 径 算 法 能 够 快速 求 得 时 变 网 络 最 短路 径 的 全 局 最 优 解 。 提 出 


的 算法 在 智能 交通 、 路 径 规划 及 无 人 车 等 领域 有 更 加 广阔 的 应 
用 前 景 。 
表 7 美国 道路 地 图 下 的 实验 结 
Dijkstra Thorup (2004) 提出 的 方法 
数据 编号 
Error% CT(s) Error% CT(s) Error% CT(s) 
D, 11.37% 475.7 11.37% 2.61 0% 2.95 
De 13.48% 615.5 13.48% 3.183 0% 3.446 
D- 20.01% 3578 20.01% 8.943 0% 9.223 
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